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We present results for the low-lying spectrum of D and Dg mesons from a lattice QCD calcu- 
lation on 2+1 flavor Clover-Wilson configurations generated by the PACS-CS collaboration. 
In particular S- and P-wave states of charmed and charmed-strange mesons are explored for 
pion masses down to 156MeV. For the heavy quark, the Fermilab method is employed. In 
addition to groimd states, some excited states are extracted using the variational method. To 
check our setup, calculations of the charmonium spectrum are also carried out. For char- 
monium, the low-lying spectrum agrees favorably with experiment. For heavy-strange and 
heavy-light systems substantial differences in comparison to experiment values remain in 
channels with nearby scattering states. 



1 Introduction 

The spectrum of charmed mesons contains states for which quark model expectations 
[1] did not hold. In particular, the charmed-strange mesons D*q(2317) and Dsi(2460), 
which in the limit of an infinitely heavy quark form a pair of mass-degenerate states with 
= 2 (where / the total angular momentum of the light quark and P is parity) turn out 
to be narrow states with masses below the DK and D*K thresholds. Due to their unantic- 
ipated properties, it has been conjectured that these states are not of a simple qq nature. 
Lattice QCD is ideally suited to calculate the properties of hadrons from first principles. 
In recent years conceptional and algorithmic improvements have enabled calculations 
with light dynamical quarks in boxes of size ?»2.5fm and with fine lattice spacings. We 
present results from a recent calculation of charmonium and charmed mesons using dy- 
namical gauge configurations with light sea quarks corresponding to pion masses as low 
as 156MeV, which is a significant improvement over previous simulations which employed 
either the quenched approximation or very heavy sea quarks. The following section sum- 
marizes the technical aspects of our simulation and the last section presents selected results. 
For more details please refer to the full published results [2]. 
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2 Technicalities 

For this study 2+1 Clover-Wilson configurations generated by the PACS-CS collabora- 
tion [3] have been used. While the mass of the strange quark is kept close to its physical 
mass, pions made from light up and down quarks range from 702 MeV to 156 MeV. The 
number of lattice points is 32'^ x 64 and the lattice spacing has been determined [3] to be 
0.907(13)fm. Table 1 shows some of the run parameters and the number of configurations 
used. 



Ensemble 








#configs D/Ds 


1 


1.52617 


0.13700 


0.13640 


200/200 


2 


1.52493 


0.13727 


0.13640 


-/200 


3 


1.52381 


0.13754 


0.13640 


200/200 


4 


1.52327 


0.13754 


0.13660 


-/200 


5 


1.52326 


0.13770 


0.13640 


200/348 


6 


1.52264 


0.13781 


0.13640 


198/198 



Table 1 : Run parameters for the PACS-CS lattices [3]. All gauge configurations have been 
generated with the inverse gauge coupling /5 = 1.90 and the light quark clover coefficient 
= 1-715. The quantity c^^J denotes the heavy quark clover coefficient. 

For the charm quarks we use the Fermilab method [4-6]. In this approach the spin-averaged 
kinetic mass of the charmed-strange mesons is tuned to assume its physical value. Once 
the tuning has been done, differences in the rest masses are to be compared to experiment. 
The resulting charm quark hopping parameter Kc is listed in Table 1. 

The spectrum results are obtained using the variational method [7,8]. For a given set of 
quantum numbers a matrix C(f),y of interpolating fields is constructed 

(1) C(0y = i:e-'="(0|O,-|n)(n|Oj|0). 

n 

On each time slice the generalized eigenvalue problem (GEVP) is solved 

(2) C{t){p^'^ = A^'\t)C{to){p^'\ 

Asymptotically only a single state contributes to each eigenvalue. For details of our quark 
sources, which contain Jacobi-smeared [9,10] and derivative sources [11,12] please refer to 
the published results [2]. 
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Figure 1 : Left panel: Mass splittings in the charmonium spectrum compared to the spin- 
averaged ground state mass (M,^^. -|- 3Mjy) / 4. Right panel: Same for charmed mesons. 



3 Results 

In this section we present a selection of results for charmed mesons and for charmonium. 
The charmonium spectrum below the DD and D*D thresholds is a good test case for our 
setup^, as it contains only well established and uncontroversial states that are all believed 
to be predominantly of a 'qq nature and as our tuning procedure uses no input from the 
charmonium spectrum. The results we obtain for the lightest sea quarks are plotted in 
the left panel of Figure 1. Some additional results are also shown in Table 2. In general 
our charmonium results agree qualitatively with the experimental spectrum. The Spin- 
dependent splittings displayed in Table 2 are expected to be sensitive to discretization 
effects and we expect that these are the dominant reason for underestimating the splittings 
compared to experiment. 



Mass difference 


Our results [MeV] 


Experiment [MeV] 


IS hyperfine 


97.8 ± 0.5 ± 1.4 


116.6 ±1.2 


IP spin-orbit 


37.5 ± 2.4 ± 0.5 


46.6 ±0.1 


IP tensor 


10.44 ±1.13 ±0.15 


16.25 ± 0.07 


2S hyperfine 


48 ± 18 ± 1 


49 ±4 



Table 2: Spin dependent mass splitting in the Charmonium spectrum. 



The right panel of Figure 1 shows the D and Dg spectrum as extracted from our simulation 
on the ensemble with the lightest sea quark mass. We obtain reasonable values for the 
hyperfine splittings and for the pairs of states corresponding to the multiplet with f = | 

^in particular with regard to the discretization of the heavy charm quark 
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Figure 2: Measured energy levels for the D*q (left panel) and Dgj (right panel) ground 
states (black circles) compared to experimental states (magenta circles). All masses are 
plotted with respect to the spin-averaged ground state. The plus signs denote the DK 
and scattering levels on the lattice (black) and in nature (red). At our lowest pion 
mass the artificially heavy scattering states are very close to the measured groimd state 
energy. 



in the heavy quark limit. While light sea quarks are important for several of the states, the 
doublets corresponding to f = ^ in the heavy quark limit show a substantial difference 
compared to experiment, which is hard to explain by discretization effects alone. For these 
states the nearby DK and D*K thresholds may play an important role. Figure 2 compares 
the measured energy levels in the lattice simulation with the experimental resonance and 
the relevant scattering threshold from experiment. It can be seen that the scattering thresh- 
old in our simulation is slightly unphysical. While the state we observe seems to coincide 
with the scattering threshold at the smallest pion mass, this is no longer the case for larger 
pion masses and the overlaps of the state with our interpolator basis (as encoded in the 
eigenvectors of the GEVP) suggests that we see the same state at all pion masses. In light 
of this, further studies should include the relevant multi-meson states in the variational 
basis, which is challenging from a computational point of view. 
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